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ABSTRACT 

Aims. We aim to determine the physical and chemical properties of the starless and protostellar cores in Orion B9, which represents 
a relatively quiescent star-forming region in Orion B. 

Methods. We observed the NH 3 (J, K) = (1, 1) and (2, 2) inversion lines and the N 2 H + (3 - 2) rotational lines, with the Effelsberg 
100-m and APEX telescopes, respectively, towards the submillimetre peak positions in Orion B9. These data are used in conjunction 
with our APEX/LABOCA 870 jim dust continuum data of the region. 

Results. The gas kinetic temperature in the cores derived from the NH3 data is between ~ 9.4 - 13.9 K. The non-thermal velocity 
dispersion is subsonic in most of the cores. The non-thermal linewidth in protostellar cores appears to increase with increasing 
bolometric luminosity. The core masses, ~ 2 - 8 Mq, are very likely drawn from the same parent distribution as the core masses in 
Orion B North. Based on the virial parameter analysis, starless cores in the region are likely to be gravitationally bound, and thus 
prestellar. Some of the cores have a lower radial velocity than the systemic velocity of the region, suggesting that they are members 
of the "low-velocity part" of Orion B. The observed core-separation distances deviate from the corresponding random-like model 
distributions. The distances between the nearest neighbours are comparable to the thermal Jeans length. The fractional abundances of 
NH3 and N2H + in the cores are ~ 1.5 - 9.8 x 10~ 8 and ~ 0.2 - 5.9 x 10~ 10 , respectively. The NH3 abundance appears to decrease 
with increasing Hi column and number densities. The NH3/N 2 H + column density ratio is larger in starless cores than in cores with 
embedded protostars. 

Conclusions. The core population in Orion B9 is comparable in physical properties to those in nearby low-mass star-forming regions. 
The Orion B9 cores also seem to resemble cores found in isolation rather than those associated with clusters. Moreover, because the 
cores may not be randomly distributed within the region (contrary to what was suggested in our Paper I), it is unclear whether the 
origin of cores could be explained by turbulent fragmentation. On the other hand, many of the core properties conform to the picture 
of dynamic core evolution. The Orion B9 region has probably been influenced by the feedback from the nearby Ori OB lb group, and 
the fragmentation of the parental cloud into cores could be caused by gravitational instability. 

Key words. Stars: formation - ISM: clouds - ISM: individual objects: Orion B - ISM: molecules - ISM - Radio lines: ISM 



1. Introduction 

The Orion B molecular cloud (LI 630) is a useful target for 
studying dense cores and the processes of star formation. The 
cloud complex lies at a relatively close distance to the Sun 
(~ 450 pc), and it contains a wide range of star-forming envi- 
ronments, such as the high-mass star-forming region NGC 2024 
and several other regions of clustered star formation (which is 
the dominant mode of star formation in the Galaxy), but also 
more quiescent areas (see, e.g., |Bally 2Q0E\ llkeda e t al. 2009; 
Buck le et al. 20 10). We recently mapped region of the central 
part of Orion B of about 36' x 27' (4.7 pc x 3.5 pc), called 
Orion B9 (Z = 206!1, b = -15°8), at 870 pm dust continuum 
using LABOCA on APEX, and discovered 12 dense submm 
cores (Miettine n et al. 200 9; hereafter Paper I). The cores were 
classified into starless and protostellar ones by using the Spitzer 



data, yielding the result that half the cores have embedded proto- 
stars). In the present paper, we examine the physical character- 
istics of the Orion B9 cores further. To determine the gas kinetic 
temperature, kinematics, and the dynamical state of the cores, 
we performed NH3(1, 1) and (2, 2), and N2H + (3 - 2) observa- 
tions towards the core submm peak positions with the Effelsberg 
100-m and APEX telescopes, respectively. In this paper, the de- 
rived temperatures are also used to recalculate the temperature- 
dependent core parameters presented in Paper I. 

This paper is organised as follows. The observations and 
data-reduction procedures are described in Sect. 2. The obser- 
vational results are presented in Sect. 3. In Sect. 4, we describe 
the analysis of the physical and chemical properties of the cores. 
Discussion of the results is given in Sect. 5, and in Sect. 6, we 
summarise our main conclusions. 



Send offprint requests to: O. Miettinen 

* This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX) under programme 084.F-9312, and ob- 
servations with the 100-m telescope of the Max-Planck-Institut fiir 
Radioastronomie (MPIfR) at Effelsberg. APEX is a collaboration be- 
tween the MPIfR, the European Southern Observatory, and the Onsala 
Space Observatory. 



2. Observations and data reduction 

2.1. Effelsberg 100-m observations 

Pointed observations of the NH3(1, 1) and (2, 2) inversion line 
emission towards the dense cores in Orion B9 were performed 
with the Effelsberg 100-m telescope of the MPIfR on 2009 
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November 23-25. The target positions listed in Table Q] were 
drawn from the submm peak positions in the LABOCA 870 jjm 
dust continuum map of the region. In Table [TJ we also show the 
source classification from Paper I. 

The 1.3 cm HEMT receiver was tuned to a frequency of 
23 708.564 MHz, lying midway between the rest frequencies of 
the NH3(1, 1) and (2, 2) lines, which are separated by about 28 
MHz. For the sky frequency conversion the local standard of rest 
(LSR) velocity was set to 9.2 km s _I . The receiver measures or- 
thogonal linear polarisations. The sum of the two channels was 
fed into the Fast Fourier Transfrom Spectrometer (FFTS) with a 
100 MHz bandwidth. In this manner, the Stokes / spectra of the 
two lines were measured simultaneously. The backend is a mod- 
ified copy of the FFTS at the APEX telescope (see Sect. 2.2). 
The 100 MHz band was split into 16 384 channels. This resulted 
in a channel separation of 6.1 kHz which corresponds to about 
77 m s~' at 23.7 GHz. The spectral resolution (equivalent noise 
bandwidth) is about 9.8 kHz (124 m s~'). One of the sources, 
SMM 6, was also observed in NH 3 (1, 1) (23 694.4955 MHz) us- 
ing a 20 MHz bandwidth to resolve the hyperfine structure in 
more detail. In this configuration, the channel separation is 1 .22 
kHz (15 ms 1 ). 

The FWHM (full width at half maximum) beamsize at the 
observed frequencies is 40" (0.09 pc at 450 pc), and the main 
beam efficiency is 77MB = 0.52. Observations were conducted 
in frequency switching mode with a frequency throw of 6 MHz. 
The single-sideband (SSB) system temperatures were about 200- 
3 10 K in the main-beam brightness temperature, Tmb, scale dur- 
ing the observations. Typical integration time was ~ 120 min per 
position, resulting in a lcr rms noise of about 38-75 mK. For fur- 
ther information of the telescope and the receiving system, see 
http : //www . rapi f r-bonn . mpg . de/r adioteleskop/. 

The NH3(1, 1) and (2, 2) emission was detected towards 
all submm cores except IRAS 05412-0105 (rms~79 mK). This 
source was also very weak (0.17 Jy beam -1 ) in our LABOCA 
870 yum map (Paper I; Fig. 1 therein). We performed test mea- 
surements of the NH 3 (3, 3) inversion line at 23 870.1292 MHz 
towards two strong NH3(1, 1) sources in our sample: SMM 4 
and IRAS 05405-01 17. No lines were detected after 30 min in- 
tegration per position (rms ~ 200 - 300 mK at the velocity res- 
olution 15 m s _1 ). 

Telescope pointing and focus were checked about every 1- 
1.5 h by continuum scans on radio quasars PKS 0420-014 and 
3C147, and the radio galaxy 3C213. The pointing was typi- 
cally accurate to within ~ 7". Absolute flux calibration was 
based upon continuum cross scans on the planetary nebula 
NGC 7027, and the radio galaxies/quasars 3C123, 3C147, and 
3C286, for which we adopted the flux densities 5.58 + 0.10, 
2.93, 1.91 + 11, and 2.46 + 05 Jy at 1.3 cm, respectively 
dOtt et al. 19941 |Peng et al. 2000) 1. We were not able to retrieve 
the zenith opacity, t z , by doing a sky-dipping measurement or 
by observing calibrators in good weather in a wide range of el- 
evations (e.g., Pill ai et al. 20061 Frieswijk et al. 2007] l. Instead, 
we estimated the average t z at 1.3 cm by fitting all the cali- 
bration measurements with the function G(6) x e~ Tz / smS , where 
G{6) is the instrumental gain-elevation curve given on the tele- 
scope www pages (see above). In this manner, we obtained an av- 
erage t z of 0.15. This value is consistent with previous estimates 
at centrimetre wavelengths (Frieswijk et al. 2007 Appendix A 
therein). The typical uncertainty in absolute flux calibration is 
~ 15% (excluding the systematic error in t z ). Note that this un- 
certainty does not affect the parameters that depend on the line 
intensity ratios (e.g, optical thickness, kinetic temperature), or 
the kinematical parameters (centroid velocity and width of the 



Table 1. Submillimetre cores in Orion B9. 



Source 


Peak position 


Class 


name 


«20oo.o [h:m:s] 


#2000.0 [° : ' : "] 




IRAS 05399-0121 


05 42 27.4 


-01 19 50 


0/1 


SMM 1 


05 42 30.5 


-01 20 45 


prestellar 


SMM 2 


05 42 32.9 


-01 25 28 


prestellar 


SMM 3 


05 42 44.4 


-01 16 03 





TD A C C\ZAC\Z Al 1 "7 

1KA5> 0341D-01 1 / 


05 43 02.7 


-01 16 21 


u 


SMM 4 


05 43 03.9 


-01 15 44 





SMM 5 


05 43 04.5 


-01 17 06 


prestellar 


SMM 6 


05 43 05.1 


-01 18 38 


prestellar 


Ori B9 N 


05 43 05.7 


-01 14 41 


prestellar 


SMM 7 


05 43 22.1 


-01 13 46 


prestellar 


IRAS 05412-0105 


05 43 46.4 


-01 04 30 





IRAS 05413-0104 


05 43 51.3 


-01 02 50 






line). The calibration uncertainty propagates to the values of the 
line excitation temperature, and NH3 column densities and frac- 
tional abundances (see Sect. 4.1). 

The spectra were reduced using the CLASS package of the 
IRAM's GILDAS software^ For a given source, the spectra were 
averaged and folded. A first order polynomial, and in one case 
(SMM 4) a third order polynomial, was applied to correct the 
baseline in the NH3(1, 1) lines. A polynomial baseline of order 
3 was subtracted from the NH?(2, 2) lines. We fitted the hyper- 
fine structure of the NH?(1, 1) line using "method NH3(1, 1)" 
of the CLASS package to derive the LSR velocity (vlsr) of the 
emission, FWHM linewidth (Av), and the line optical thickness 
(see Sect. 3.1). The emission of the (2, 2) satellite lines was not 
detected. Nevertheless, the radial velocities and linewidths were 
determined using "method NH3(2,2)" of CLASS. 



2.2. APEX observation 

The N2H + (3 - 2) observations towards nine cores in the central 
region of Orion B9 (i.e., IRAS 05399-0121 and 05405-0117, 
SMM 1, and 3-7, and Ori B9 N) were carried out on 2009 
September 2 and 1 1 with the APEX telescop^B The backend was 
the MPIfR FFTS dKlein et al. 20061 1 with a 1 GHz bandwidth di- 
vided into 8192 channels. The resulting channel width was 122 
kHz which corresponds to 0. 13 km s _1 at the observed frequency 
279.5 GHz. As frontend, we used APEX-2 of the Swe dish 
Heterodyne Facility Instrument (SHFI; Vas silev et al. 20081 . At 
279.5 GHz, the APEX beamsize is ^ 22'.'3, and the main beam 
efficiency is t]mb — 0.74. 

The observations were performed in the wobbler-switching 
mode with a 100" throw (symmetric offsets) and a chopping rate 
of 0.5 Hz. The typical total integration time was ~ 1 1 min per 
position, and the SSB system temperature was between 210-360 
K (Tmb scale) during these observations. The telescope pointing 
was checked using the planets Mars and Uranus, and the stars a 
Orionis and V1259 Ori, and was found to be accurate to ~ 3". 
The calibration was achieved by the chopper-wheel method, and 
the intensity scale given by the system is T* A , the antenna temper- 
ature corrected for atmospheric attenuation. The observed inten- 
sities were converted to the main-beam brightness temperature 
scale by T MB = T*Jj] MB . 

The spectra were reduced using the CLASS. The individual 
spectra were averaged, and linear baselines were subtracted from 



1 http://www.iram.fr/IRAMFR/GILDAS 

2 Other spectral lines observed under this programme include 
C n O(2 - 1), N 2 D + (3 - 2), DCO + (4 - 3), and H 13 CO + (4 - 3). These 
observations will be presented in the forthcoming paper. 
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the resulting sum spectra. The resulting lcr rms noise values are 
35-130 mK. 

The 7 = 3-2 transition of N2H + contains 38 hyperfine 
components. We fitted the hyperfine structure of the lines using 
"method hfs" of the CLASS. For the rest frequencies of the hy- 
perfine components, we used the values from Pagani et al. (2009; 
Table 4 therein). The adopted central frequency, 279 511.832 
MHz, is that of the Jf t f = 345 — > 234 hyperfine component 
which has a relative intensity of 17.46 % 

3. Observational results 

3.1. NH 3 

The NH3(1, 1) and (2, 2) spectra are shown in Fig. Q] The 
NH 3 (1, 1) spectrum observed from SMM 6 with a 20 MHz 
bandwidth (see Sect. 2.1) is shown in Fig. [2] In Table [2] we 
give the line parameters resulting from the hyperfine/Gaussian 
fits to the NH3 lines. The so-called 'main group optical thick- 
ness', r m , given in Table|2]is the sum of peak optical thicknesses 
of the eight hyperfine components included in the main group. 
For the (1,1) transition this equals to the corresponding sum 
of the ten satellite components, so the total of all peak opti- 
cal thicknesses is 2 x r m . Assuming that the line profile (due 
to velocity dispersion) is Gaussian, the integrated optical thick- 
ness can be obtained from J r(v)dv = V"V hi2r m Av, where 
Av is the linewidth (FWHM) in velocity units. The uncertain- 
ties associated with Vlsr, Av, and r m are the lcr fitting errors. 
The uncertainty in 7mb, 0"(7mb)> includes the absolute calibra- 
tion uncertainty of 15%, and the lcr rms noise in the spectrum, 
and was obtained by quadratic summing of these two errors as 

c(7mb) = -^cr^ al + cr;? ms . The determination of the line excita- 
tion temperature, T ex , listed in the last column of Tabled is de- 
scribed in Sect. 4.1.1. 

The NH 3 spectra towards SMM 4 and Ori B9 N show evi- 
dence of a second velocity component in addition to the strong 
line at about 9.1 km s _1 corresponding to the systemic velocity 
of Orion B. In the case of IRAS 05405-01 17, there is also a hint 
of a second velocity component, probably resulting from the fact 
that the 40" beam slightly overlaps with the nearby core SMM 

4. In the direction of SMM 4, the second component has a veloc- 
ity of about 1.6 km s _1 , whereas towards Ori B9 N the velocity 
is 1.9 km s _1 . In both cases, the main hyperfine group of the 
second velocity component overlaps with the inner satellite line 
F\ = 1 — 2 of the principal velocity component. Consequently, 
the inner satellite line F\ — 2 - 1 of the second velocity com- 
ponent contaminates the main hyperfine group of the principal 
velocity component. These spectra were analysed using a two- 
component fit to the NH3(1, 1) hyperfine structure. We note that 
the secondary velocity components towards SMM 4 and Ori B9 
N are close to the radial velocities of IRAS 05413-0104 (1.5 
km s" 1 ) and SMM 7 (3.6 km s ). We also note that the second 
velocity components influence the core parameters derived from 
the dust continuum emission (see Sect. 4.3). Part of the observed 
flux density may be due to the dust component associated with 
another source along the line of sight. However, because con- 
tinuum observations do not provide velocity information, it is 
impossible to solve this "overlap" problem. 

The NH3(1, 1) line towards SMM 6 measured using a 20 
MHz bandwidth is 30 m s narrower and has a 15% larger op- 
tical thickness than that measured using a 100 MHz bandwidth. 
The difference in the linewidths corresponds to effect expected 
from instrumental broadening with the two configurations. As is 



expected, the product Avr m is roughly constant at the two spec- 
tral resolutions (see Eq. ©). As the high resolution spectrum is 
only available for one object, we use in the subsequent analysis 
the 100 MHz spectra for all of them. 



- SMM 6 

\ NH 3 (1 ,1) (20 MHz aw) 
- 

Fj = l-0 F i = 1 - 2 

1 1 1 


U 

1 1 


1 1 

uJ 

= 2-1 F, 

1 1 


1 obs. " 
=0-1 


residuals (obs. — fit) 

I , , , I 



Fig.2. NH3(1, 1) spectrum measured from SMM 6 using the 
FFTS with a 20 MHz bandwidth. Below the spectrum are shown 
the model fit to the hyperfine structure and the residual spec- 
trum. Satellite hyperfine components in the Fj — > F\ format are 
indicated. 



3.2. N 2 H + 

The N2H + (3 - 2) spectra are shown in Fig. [3] In Table [3] we 
show the line parameters obtained from the fits to the hyperfine 
structure. Here, the value of <t(Tmb) includes only the lcr rms 
noise in the spectrum (Col. (4) of Table[3]). 

The LSR velocities of the N2H + (3 - 2) lines are mostly simi- 
lar to those determined from NH3(1, 1), as can be seen in the top 
panel of Fig. [4] On the other hand, there is hardly any correla- 
tion between the N 2 H + (3 - 2) and NH 3 (1, 1) linewidths (bottom 
panel of Fig . |4j» . It should be noted, however, that the hyperfine 
fitting to the N 2 H + spectra is uncertain due to strongly overlap- 
ping components and rather a poor signal-to-noise (S/N) ratio. 

The N 2 H + (3 - 2) spectra towards SMM 4 and Ori B9 N show 
the same lower radial velocity components (< 2 km s _1 ) as seen 
in NH3(1, 1). On the other hand, the 'principal component' at 
~ 9 km s^ 1 is not detected towards SMM 4 and it is also very 
weak towards Ori B9 N. Similarly, in Paper I we found that the 
N 2 H + (1 - 0) spectra towards the selected position near IRAS 
05405-01 17/SMM4 and Ori B9 N have second velocity compo- 
nents at ~ 1.3 and ~ 2.2 km s _1 , respectively. The latter position 
also had additional N 2 D + (2 - 1) line centred at ~ 2.3 km s . 

We note that Harju et al. (2006) probably detected N 2 H + (4 - 
3) towards IRAS 05405-0117. They were using a double side- 
band system where a coincidence with DCO + (5 - 4) from im- 
age band could not be ruled out. Using the rest frequencies of 
N 2 H + rotational lines from Pagani et al. (2009), the peak veloc- 
ity of the suggested N 2 H + (4 - 3) line (9.25 km s _1 ) is very sim- 
ilar to that of N 2 H + (3 - 2) (9.30 km s 1 ). For DCO + (5 - 4) the 
peak velocity would be about 0.3 km s lower (8.95 km s ). 
The N 2 H + (4 - 3) linewidth (0.34 km s _1 ) is similar to that of 
NH 3 (1, 1) (0.32 km s- 1 ). 
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IRAS 06,595—01 ^ ' 



NHj(2,2) ft 



;NH 3 (2,2) . : 



NH 3 (1,1) 



- 1 ^f^yXyJh^ 



40 -20 



; I I || | 



RAS 05405-01 17 



40 -20 



NH 3 (2,2) 




SMM 6 

NH 3 (2,2) 








40 -20 



NH 3 (2,2) , 



NHj(2,2) 



RAS (!!.Vl'2-G105 



^NH 3 (2,2; 



HAS 05413-0104 




A NH 3 {2,2) 



Fig. 1. NH3(1, 1) and (2, 2) spectra measured from the pre- and protostellar cores in Orion B9. The temperature scale is in Tmb- 
NH3(2, 2) spectra are offset by -1 K from zero baseline for clarity. Note that there are two velocity components towards SMM 4 
and Ori B9 N, and that no lines were detected towards IRAS 05412-0105. 



4. Analysis 

4.1. NH 3 analysis 

In this subsection, we derive physical parameters from NH3 data. 
The obtained results are given in Col. (7) of Table[2](r ex ), and in 
Table [4] The analysis follows the concept first presented in Ho 
et al. (1979) and further discussed in Ho & Townes (1983). Here 
it is assumed that the same excitation temperature, r ex , charac- 
terises all the hyperfine transitions between the (1, 1) and (2, 2) 
inversion doublets. 

4.1.1. Excitation, rotation, and kinetic temperatures 

The T ex was calculated at the (1,1) line peak using the optical 
thickness, T pe ak, at that velocity, from the equation 



(1) 



ln[l 



7"mb(v) 1 



F(T bg )] 



where T\\ = Iivu/Icb, h is the Planck constant, vn is the (1, 1) 
line frequency, kg is the Boltzmann constant, 7b g = 2.73 K is 

the background temperature, and F(T) — (e, hv l k * T - l) . The 
uncertainty in T ex was calculated by propagating the errors in 
the peak values of 7mb(v) and r(v). The r(v) distribution was 



obtained using the (1, 1) hyperfine fit from CLASS which yields 
the LSR velocity line centroid, the width of an individual hyper- 
fine component, and the main group optical thickness, T m (see 
Sect. 3.1). 

The rotational temperature, T rot , can be in principle calcu- 
lated using the formula 



-41.5 



m h A<(i,i)J 



(2) 



where N(l, 1) and N(2, 2) are the column densities of ammo- 
nia molecules at the (1, 1) and (2, 2) rotational levels, respec- 
tively. The estimation of these column densities is described 
in Sect. 4.1.2. The ratio 3/5 refers to the ratio of the statistical 
weights of the two states (g(l, l)/g(2, 2)). 

When the optical thickness, r m , is known, Eq. (f2]) can be 
written as 



Trnt — 



-41.5 



In 



-0.283 



ln[l 



r MB (2.2) 



(1 



■)]} 



(3) 



Jmb(1, 1) 

dHo & Townes 19831 their Eq. (4)). Here it is assumed that T ex 
and Av are the same for both the (1, l)and(2, 2) transitions. The 
error associated with T ml was propagated from the uncertainties 
in TmbO, 1), 7mb(2, 2), and r m . When T m could not be reliably 
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Table 2. NH 3 (1, 1) and (2, 2) line parameters. 





Transition 


VLSR 


Av 


7 a 


T h 


T 


Source 


(J,K) 


[km s _1 ] 


[kms" 1 ] 


[K] 




[K] 



IRAS 05399-0121 


(1. 


1) 


8.61 


±0.01 


0.59 


±0.02 


1.6 


±0.2 


1.10 


±0.13 


5.4 ± 


0.5 




(2, 


2) 


8.64 


±0.02 


0.42 


±0.08 


0.3 


±0.1 










SMM 1 


(1. 


1) 


9.04: 


t 0.003 


0.67 


±0.01 


2.9 


±0.4 


1.23 


±0.06 


7.2 ± 


0.7 




(2. 


2) 


9.13 


±0.02 


0.53 


±0.05 


0.4 


±0.1 










SMM2 


(1. 


1) 


9.71 : 


t 0.003 


0.34 


±0.01 


1.4 


±0.2 


1.38 


±0.16 


5.3 ± 


0.4 




(2, 


2) 


9.72 


±0.02 


0.26 


±0.04 


0.2 


±0.1 










SMM 3 


(1, 


1) 


8.40: 


t 0.002 


0.38 


±0.01 


2.5 


±0.4 


2.14 


±0.10 


6.1 ± 


0.5 




(2, 


2) 


8.42 


±0.01 


0.42 


±0.06 


0.4 


±0.1 










IRAS 05405-0117 


(1. 


1) 


9.12: 


t 0.001 


0.32: 


t 0.003 


4.9 


±0.7 


3.36 


±0.09 


8.4 ± 


0.9 




(2. 


2) 


9.17 


±0.01 


0.30 


± 0.02 


1.1 


±0.2 










SMM 4 


(1, 


1) 


9.12: 


t 0.001 


0.34 


±0.01 


3.2: 


±0.5^ 


1.33 


±0.03 


8.6 ± 


0.9 




(2, 


2) 


9.15 


±0.01 


0.23 


±0.02 


0.7 


±0.f 










SMM 4 (2nd v-comp.) 


(1, 


1) 


1.62: 


t 0.003 


0.71 


±0.01 


2.0 


±0.2 


1.50 


±0.06 


5.1 ± 


0.5 




(2, 


2) 


1.66 


±0.03 


0.51 


± 0.07 


0.2 


±0.1 










SMM 5 


(1. 


1) 


9.16: 


t 0.002 


0.35 : 


t 0.004 


3.0 


±0.5 


2.28 


±0.08 


6.6 ± 


0.6 




(2. 


2) 


9.21 


±0.01 


0.31 


± 0.03 


0.5 


±0.1 










SMM 6 


(1. 


1) 


9.38: 


t 0.001 


0.29: 


t 0.002 


3.4 


±0.5 


4.66 


±0.11 


6.3 ± 


0.5 




(2, 


2) 


9.41 : 


t 0.004 


0.30 


±0.01 


0.9 


±0.1 










SMM 6 (20 MHz BW) 


(1. 


1) 


9.44: 


t 0.001 


0.26: 


t 0.002 


3.3 


±0.5 


5.38 


±0.13 


6.2 ± 


0.5 


Ori B9 N 


(1. 


1) 


9.14: 


t 0.004 


0.40 


±0.01 


1.1 : 


±0.2 C 


1.08 


±0.15 


5.0 ± 


0.4 




(2, 


2) 


9.18 


±0.02 


0.21 


±0.04 


0.2: 


t0.05 










Ori B9 N (2nd v-comp.) 


(1. 


1) 


1.93 


±0.01 


1.50 


±0.04 


0.7 


±0.1 


0.37 


±0.09 


5.1 ± 


0.6 




(2, 


2) 


2.29 


±0.11 


2.53 


±0.28 


0.1 : 


t0.04 










SMM 7 


(1. 


1) 


3.60 


±0.01 


0.60 


±0.02 


1.2 


±0.2 


2.14 


±0.14 


4.1 ± 


0.2 




(2. 


2) 


3.89 


±0.10 


0.67 


±0.14 


0.1 : 


t0.05 










IRAS 05413-0104 


(1. 


1) 


1.48: 


t 0.003 


0.46 


±0.01 


1.8 


±0.3 


1.68 


±0.09 


5.2 ± 


0.4 




(2, 


2) 


1.52 


±0.02 


0.53 


±0.04 


0.4 


±0.1 











" The error in Tmb includes the ~ 15% calibration uncertainty (see the text in Sect. 3.1). 

b T m is the optical thickness of the (1, 1) main group. The total optical thickness, i.e., the sum over all satellites, is T tc 
c 7mb(1, 1) was determined by substracting the second velocity component from the spectrum. 



2xt, t 



Table 3. N 2 H + (3 - 2) line parameters. 



Source 




V LSR 

[km s" 1 


] 


Av 
[kms" 1 ] 




r a 
1 

[K] 


T 


. b 

tot_ 


T 

' ex 

[K] 




IRAS 05399-0121 


8 


.73 — 


= 


.02 


0.67 = 


t0.07 


0. 


.92 


±0.04 


3.08 


± 1.10 


5.1 ± 0. 


3 


SMM 1 


9 


.30 = 


= 


,03 


0.61 = 


tO.09 


0. 


.66 


±0.07 


2.77 


± 1.48 


4.6 ± 0. 


3 


SMM 3 


8 


.57 = 


= 0. 


.03 


0.71 = 


= 0.02 





.65 


±0.08 


0.92 


±0.19 


5.5 ±0. 


5 


IRAS 05405-0117 


9 


,30 i 


= 


,03 


0.68 = 


tO.06 


0. 


66 


±0.07 


1.48 


± 0.34^ 


5.q£ 




SMM 4^ 




























SMM 4 (2nd v~comp.) 


1 


.60 = 


= 0. 


.10 


0.44 = 


= 0.25 


0. 


.35 


±0.13 


0.54 


± 0.26 1 ' 


5.0 rf 




SMM 5 


9 


,43 = 


= 0. 


.07 


0.46 = 


= 0.13 


0. 


.17 


±0.06 


0.23 


± 0.09 c 


5.0 rf 




SMM 6 


9 


.52 = 


= 


.03 


0.58 = 


tO.09 


0. 


.53 


±0.08 


0.99 


± 0.24 c 


5.0 rf 




OriB9N 


8 


,78 = 


= 0. 


.13 


0.22=1 


: 0.32£ 




< ' 


0.07 


< 0.10 1 ' 


5.0 rf 




Ori B9 N (2nd v-comp.) 


2 


,07 = 


= 


.10 


0.43 = 


t0.22~ 


0. 


.10 


±0.05 


0.13 


± 0.07 c 


5.0 rf 




SMM 7 


3 


.96 = 


= 0. 


.14 


0.66 = 


tO.22 





.13 


±0.08 


0.17 


±0.11 c 


5.0 rf 





" The error in Tub is the lcr rms noise error. 

* T lot is the total optical thickness of the hyperfine multiplet. 

c T( 0t could not be determined through fitting the hyperfine structure. Instead, T tol was calculated by assuming T a 
account that the line intensity in Col. (4) accounts for 92.6% of the total line strength. 
d Assumed value. 

' The principal velocity component was not detected. 
' The associated error is larger than the value. 



5.0 K, and taking into 



determined (SMM 4, Ori B9 N, and their additional velocity 
components), we estimated the NH3(1, 1) and NH3(2, 2) col- 
umn densities from the integrated intensities assuming T ex = 5 
K, as described below. In these cases, r iot was also calculated 
using Eq. (0, and the uncertainty in T TOt was propagated from 
the uncertainties in the column densities. The values of T Iot cal- 



culated using Eq. are mostly similar to those resulting from 
Eq. (01 within the errors (see Table |4}. The subsequent analy- 
sis includes only those T Iot values which could be derived using 
Eq. ©. 

The gas kinetic temperature, 7\j n , was calculated from T rot 
using the relationship given by Tafalla et al. (2004): 




Fig. 3. N2H + (3 - 2) spectra measured from the pre- and protostellar cores in Orion B9. Spectra are overlaid with the 38 component 
hyperfine structure fits. Hyperfine fits to the second velocity components are indicated by red lines. Note that only the second 
velocity component is detected towards SMM 4. The relative velocity of each individual hyperfine component is labelled with a 
short bar on the spectrum towards IRAS 05399-0121 (top left panel). The temperature scale is in Tmb- 



7kii 



1 



(4) 



1 + Lie r.™ 



This relationship is recommended for dense cores with gas tem- 
peratures < 20 K. Uncertainty in T^ a was propagated from the 
uncertainty in T rot . 

4.1.2. NH 3 column density calculations 

The column density in the (J, K) — (1, 1) state was calculated 
using the formula (see Harju et al. 1993 Eqs. (3) and (5) therein) 



3/t «o 



^ y(7+ ^(r ex )(e^^ + l)AvT t( 



2^ 2 VhT2 K 
7.83 x 10 n F(T ex ) (e U4/r « + l) Av[km s^]^, cm 



(5) 



where eo is the vacuum permittivity, fi is the permanent electric 
dipole moment (1.476 D), and r tot = 2 x r m . 

For SMM 4 and Ori B9 N, and their second velocity com- 
ponents (for which r m is uncertain in all the cases), N{1, 1) was 
also calculated by assuming optically thin emission (r «: 1), 
and using the integrated intensities of clean satellites, i.e., those 
where the two velocity components do not overlap. There are 
three satellite lines in the NH 3 (1, 1) spectra of both SMM 4 and 
Ori B9 N that do not suffer from the contamination of the second 
velocity component, i.e., F\ = 1 -0, 2 - 1, and - 1 (see Fig. [I). 
For the second velocity components, the corresponding satellite 
lines are F\ = 1 - 0, 1 - 2, and - 1 . Taking into account that 
these lines comprise R, = 36.1% of the total line strength, and 



assuming that the filling fraction of the emission in the beam is 
rjf = 1 , we can derive the formula 



N(l, 1) = 



3k B e J{J+ l)e r " / 7 ° + 1 1 

F(T CX ) 



2n 2 fi 2 v u K 2 i _ £Ek2 R 



T MB {1, l;s)dv 



= 6.47 x 10 



12e M4/r "+l 1 



1- 



jj: j 7mb(1, l;s)dv[Kkms- 1 ]cm- 2 , (6) 



where Tmb is integrated over the three above mentioned satel- 
lite lines (the integrated satellite intensities are 2.53 ± 0.03 and 
0.74 ± 0.03 Kkm s~' for SMM 4 and Ori B9 N, respectively; for 
the corresponding second velocity components, the values are 
2.18 ± 0.02 and 1.18 ± 0.03 K km s" 1 , respectively). For this cal- 
culation, the excitation temperature was assumed to be 5 K for 
both SMM 4 and Ori B9 N (and their second velocity compo- 
nents). This assumption is similar to the values calculated from 
Eq. ([1} in all the other cases except the principal velocity com- 
ponent towards SMM 4, where T ex might be > 5 K (see Col. (7) 
of Table 0. Similarly, N{2, 2) (needed in the calculation of r rot 
using Eq. (f2]l) was calculated by assuming that r <§c 1, and us- 
ing the integrated intensity of the main (2, 2) hyperfine complex, 
which accounts for/?, = 79.63% of the total line strength. In this 
case, we can write 



N(2, 2) = 



3k B eo J(J+l)e 



27tVv 2 2 K 2 



• T *' T ~ + l l r 

i _ r J 

F(7"cx) 



T MB (2, 2; m)dv 



4.85 x 10 



12e 114/r "+l 1 



F(T b „) 



£ / r MB (2, 2; m)dv[K km s" 1 ] cnr 2 ,(7) 
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NH3 column densities derived for each source from Eqs. © and 
(0 to keep the data set homogeneous (see Col. (4) of TableHJ). 



4.1.3. Fractional NH 3 abundance 

The fractional NH3 abundance was calculated by dividing the to- 
tal NH3 column density by the H2 column density, i.e., x(NH3) = 
Af(NH3)/Af(H2). The H2 column densities were determined from 
the submm dust continuum emission (Paper I; Eq. (3) therein), 
using the kinetic temperatures derived in the present paper and 
assuming that 7\j n = Jdust (see also Sect. 4.3). We smoothed the 
LABOCA 870 jum map to correspond the 40" resolution of the 
Effelsberg NH3 observations. 



v lsr (NH 3 ) [km s" 1 ] 



4.1.4. Non-thermal velocity dispersion and the level of 
internal turbulence 



0.0 : 

<n -0.2 j 

E : 
* -0.4 : 

^ -0.6 
r 

S -0.8 : 

> 

^ -1.0 - 



o 



-1 .2 



-0.6 



-0.5 -0.4 
log[Av(NH, 



-0.3 -0.2 
[km s ~ | 



Fig. 4. Top: The centroid velocity of the N2H + (3 - 2) emission 
versus the velocity of the NH3(1, 1) emission. Bottom: Log-log 
plot of the N 2 H + (3-2) linewidth versus the NH 3 (1, 1) linewidth. 
Prestellar cores are indicated by open diamonds and protostellar 
cores are marked by filled diamonds. In both plots, the dashed 
line shows the equality. 



where V22 is the (2, 2) line frequency, and T22 = hv^ilka- The 
integrated intensities of the (2, 2) main line for SMM 4 and Ori 
B 9 N are 0.21 + 0.01 and 0.05 + 0.01 K km s _1 , respectively 
(0.15 ± 0.01 and 0.21 ± 0.03 K km s _1 for the corresponding 
second velocity components). The uncertainties associated with 
N{\, 1) and N(2, 2) were propagated from the formal errors in 
the integrated intensities. 

The total NH3 column density, iV(NH3), was calculated 
by scaling the column density in the (1,1) state by the 
ratio N(NHs)/N(J, 1) obtained from the partition function 



Z. This yields approximately (e.g., Rosolowsky et al. 2008 
Busquet et al. 2009[) 



AT(NH 3 ) = 7Y(1, 



l)|-e*« 



3 



41.5 

Trot + 



14 



(8) 



The uncertainty in A^NIT;) was calculated by propagating the er- 
rors associated with T ex , Av, T m , and r iot . The values of A^NHj) 
calculated by substituting N(l, 1) from Eq. (© into Eq. ([8]) are 
mostly in good agreement, within the errors, with those result- 
ing from Eqs. © and ([8]). In the subsequent analysis, we only use 



The measured NH3(1, 1) linewidths were used to calculate the 
non-thermal portion of the line-of-sight velocity dispersion (av- 
eraged over a ~ 40" beam), and the level of internal turbu- 
lence. The observed velocity dispersion is related to the FWHM 
linewidth as <x b s = Av/ V8 In 2. The non-thermal velocity dis- 
persion can then be calculated as follows: 



cr NX 



(9) 



where /j. mo \ is the mass of the emitting molecule in units of 
atomic mass number (17 for NH3), and ;«h is the mass of the 
hydrogen atom. Furthermore, the level of internal turbulence is 
given by / tul -b = <x NT /c s , where c s is the one-dimensional isother- 
mal sound speed (0.19 km s _1 in a 10 K H2 gas with 10% He). 
The errors in cr NT and / tu ,b were derived by propagating the er- 
rors in Av and 7\i n . 



4.2. N 2 H + analysis 

In this subsection, we derive physical parameters from N2I-F 
data. The obtained results are given in Col. (6) of Table[3](r ex ), 
and in Table [5] 



4.2.1 . T" ex , <t nt , and / mr b 

The excitation temperature of the N2H + (3 - 2) transition was 
calculated as in the case of NITjfl, 1) (see Eq. ([1)). The opti- 
cal thickness, and thus T ex , could be determined only for three 
sources: IRAS 05399-0121, and SMM 1 and 3. For the rest of 
the sources, it was assumed that T ex = 5 K, and T tot was then 
estimated using this value. The value T ex = 5 K is expected to 
be a reasonable choice because for the above three sources r ex 
is around 5 kelvins. 

The values of <x NX and / tur b were calculated as in the case of 
NH3(1, 1) (see Eq. (O), by using the values of T^ n from NH3 
measurements. For N2FF, yU mo i is 29. 



4.2.2. N 2 H + column density and fractional abundance 

The N2H + column density was calculated using the formula (see, 
e.g., Paper I) 
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Table 4. Parameters derived from NHU data. 





T m 




AT(NH 3 ) 


x(NH 3 ) 




f{mb 


Source 


[K] 


[K] 


[10 14 cm- 2 ] 


[10- 8 ] 


[kms- 1 ] 





IRAS 05399-0121 


12.5 


± 1.4 


13.5 


± 1.6 


3.2: 


t0.6 


1 


5 i 


:0 


4 


0.24 


±0.01 


1.1 


±0.1 


SMM 1 


11.2 


±0.8 


11.9 


±0.9 


6.2: 


h0.9 


3 


5 : 


:0 


7 


0.27 4 


: 0.004 


1.3 


±0.1 


SMM2 


11.1 


± 1.6 


11.8 


± 1.8 


2.6: 


t0.7 


4 


1 = 


zO 


7 


0.12 


±0.01 


0.6 


±0.1 


SMM 3 


10.7 


±0.7 


11.3 


±0.8 


5.6: 


h0.8 


1 


9 = 


zO 


4 


0.14 4 


: 0.003 


0.7 4 


:0.04 


IRAS 05405-0117 


10.7 


±0.6 


11.3 


±0.6 


10.0 


± 1.3 


9 


8; 


z 1 


6 


0.11 4 


: 0.002 


0.6 4 


:0.02 


SMM 4 


12.8 


±0.6 


13.9 


±0.8 


3.4: 


t0.4 


3 


8d 


:0 


6 


0.124 


: 0.005 


0.5 4 


:0.03 




10.2 


±0.9,1 


10.7 


± 1.1" 


8.8± 


1.4* 


















SMM 4 (2nd v-comp.) 


10.0 


± 1.2 


10.4 


± 1.4 


6.8: 


t 1." 


4 


5 d 


= 1 


6 


0.29 


±0.01 


1.5 


±0.1 




9.7 4 


: 1.0° 


10.2 


t 1.1" 


8.2 ± 


1.6* 


















SMM 5 


10.7 


±0.6 


11.3 


±0.7 


5.9: 


t0.7 


8 


Od 


= 1 


3 


0.13 4 


: 0.002 


0.6 4 


:0.02 


SMM 6 


10.4 


±0.3 


11.0 


±0.4 


9.9: 


h0.9 


7 


2d 


:0 


8 


0.10 4 


: 0.001 


0.5 4 


:0.01 


Ori B9 N 


12.4 


± 1.1 


13.4 


± 1.3 


2.0: 


t0.4 


3 


8 = 


zO 


9 


0.15 


±0.01 


0.7 4 


:0.04 




9.7 4 


:3.6" 


10.2 


±4.0" 


2.8± 


1.9* 


















Ori B9 N (2nd v-comp.) 


12.5 


±2.1 


13.6 


±2.5 


2.5: 


t0.8 


5 


Od 


z 2 


3 


0.63 


±0.02 


2.9 


±0.3 




12.5 


± 1.7" 


13.6 


± 1.9" 


3.0± 


0.5* 


















SMM 7 


9.0: 


t 1.0 


9.4: 


t 1.1 


8.0: 


t2.0 


5 


3d 


= 1 


9 


0.25 


±0.01 


1.3 


±0.1 


IRAS 05413-0104 


12.3 


± 1.0 


13.4 


± 1.2 


3.7: 


t0.5 


3 


4d 


:0 


7 


0.18 4 


: 0.004 


0.8 4 


:0.04 



Calculated by using the Eq. ©. 
Calculated by using Eqs. (|6j and {8j. 



N(N 2 m = 3/zeo ^ l E j kB T« F{Tex)ZiTes)AvTtot 
2n 2 fi 2 2 Vrn2 Ju 
= 2.46 x 10 11 e 26 83/7 '° J F(r ex )Z(r ex )Av[km s -1 ]^ cm" 2 , (10) 

where fj. is 3.4 D, the upper rotational level number J u = 3, E a = 
hBJ u (J u + 1) is the energy of the upper transition state, and B = 
46 586.8713 MHz is the rotational constant (Pa ganietal. 2009| ). 
The rotational partition function was approximated by Z(T ex ) = 

\ hB ^ 3J- 

In those cases where the line was optically thin (i.e., SMM 5, 
Ori B9 N and its second velocity component, and SMM 7), the 
N2H + column density was also calculated from the integrated in- 
tensity. The integrated intensities obtained from Gaussian fits are 
0.15 + 0.03, 0.07 + 0.16 (0.12 + 0.04), and 0.22 + 0.07 K km s 1 
for SMM 5, Ori B9 N (and its second velocity component), and 
SMM 7, respectively. Note that the error is larger than the value 
for Ori B9 N. By combining Eqs. (4) and (6) of Paper I, we get 
the following formula for the A^(N2H + ) as a function of inte- 
grated intensity: 



^ (N2H+) = ^J^'^^V^l I rMBdv 

= 1.72 x 10 10 Z(r ex )S^ f r MB dv[Kkm s" 1 ] cirT 2 . (11) 

When applying Eq. (TTTT i. it was taken into account that the de- 
tectable emission feature with 24 hyperfine components contains 
92.6% of the total line strength. The values of N(N 2 H + ) calcu- 
lated from Eqs. ( TTUb and (fTTb are similar to each other within the 
errors (see Table|5]). 

The fractional N2H + abundance was calculated in a simi- 
lar manner as for NH3. In this case, the LABOCA map was 
smoothed to the resolution of the N2H + observations (22'.'3). 



4.3. Revision of core properties presented in Paper I 

In Paper I, we assumed that the dust temperature is Tjust = 10 
K for the starless cores, whereas for protostellar cores we used 
the temperatures derived from the spectral energy distribution 
(SED) fits (Table 6 in Paper I). By using the gas kinetic temper- 
atures determined in the present paper, we recalculated several 
parameters presented in Paper I, by assuming that T^\ n = 7d us t- 
These include the core mass, M, H2 column density, A^Fb), and 
the average H2 number density, (n(H2)). In the present paper, we 
have also amended some formula used in the derivation of the 
physical parameters: in the calculation of (n(H2)), we have used 
the effective radius /? e ff = ylA/ji (instead of the FWHM radius 
used in Paper I), and the value 2.8»2h for the mean molecular 
weight per H2 molecule (instead of 2.33hjh which is the mean 
particle weight in an H2 + 10% He mixture). The former change 
is motivated by the fact that the flux densities used to calculate 
the core masses refer to their projected areas, A. The results of 
these calculations are presented in Table [6] The revised values 
of M and A^(H2), are, on average, about 11% and 9% higher 
than those reported in Paper I. On the other hand, the new av- 
erage densities, (n(H2)), are only about 1/5 of those calculated 
in Paper I. Note that the total mass of the cores is still ~ 50 M©, 
i.e., about 4% of the total mass in the region (Paper I). 

In Paper I, we compared the core mass distributions in Orion 
B9 and Orion B North studied by Nutter & Ward-Thompson 
(2007) (Sect. 5.2 and Fig. 7 in Paper I). The core mass func- 
tions (CMFs) were constructed by removing the Class I proto- 
stellar cores from the samples (i.e., IRAS 05399-0121 was re- 
moved from the Orion B9 core sample). We performed a two- 
sample Kolmogorov-Smirnov (K-S) test between the two CMFs 
to examine if they represent the subsamples of the same un- 
derlying parent distribution. For this purpose, the masses were 
scaled to compensate for the different assumptions about the 
distance, dust temperature, and opacity. We found a very high 
likelihood of ~ 95% for the null hypothesis that the two CMFs 
are drawn from the same parent distribution (i.e., the probabil- 
ity that the two samples do not have the same parent distribution 
is 1 - 0.95 = 0.05). We repeated this analysis by using the up- 
dated masses listed in Col. (3) of Table|6] In this case, a K-S test 
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Table 5. Parameters derived from N2H + (3 - 2) data. 





N(N 2 H + ) 


x(N 2 H + ) 


N(NH 3 )/N(N 2 ll + ) 


tr NT 


/turb 


Source 


[10 13 cm- 2 ] 


[10~ 10 ] 




[kms 4 ] 






IRAS 05399-0121 


2.0 ±0.7 


4.3 ± 1.8 


16 ±7 


0.28 ± 0.03 


1.3 ± 


0.2 


SMM 1 


1.9 ± 1.1 


5.9 ±3.4 


32 ± 18 


0.25 ± 0.04 


1.2 ± 


0.2 


SMM3 


0.6 ±0.1 


0.8 ±0.2 


99 ±25 


0.30 ±0.01 


1.5 ± 


0.1 


IRAS 05405-0117 


1.0 ±0.2^ 


5.5 ± 1.5 


100 ± 28 


0.28 ± 0.03 


1.4 ± 


0.1 


SMM 4^ 


- 


- 


- 


- 


- 




SMM 4 (2nd v~comp.) 


0.2 ±0.2 


0.8 ±0.6 


288 ± 225 


0.18 ±0.11 


0.9 ± 


0.6 


SMM 5 


0.1 ±0.1 


0.7 ± 0.4 


562 ± 279 


0.19 ±0.06 


0.9 ± 


0.3 




0.2 ± 0.03^, 


1.2 ± 0.2 C 


275 ± 54' 








SMM 6 


0.6 ± 0.2~ 


2.1 ±0.6 


194 ± 80 


0.24 ± 0.04 


1.2 ± 


0.2 


Ori B9 N 


0.02 ± 0.03 


0.2 ± 0.3 


916 ± 1345^ 


0.07 ±0.18'' 


0.3 ± 


0.8 




0.1 ±0.2 c - rf 


0.8±0.1 c 


280 ± 190~ 








Ori B9 N (2nd v-comp.) 


0.1 ±0.04 


0.6 ±0.5 


250 ± 90 


0.17 ±0.10 


0.8 ± 


0.5 




0.1 ±0.05' 


1.1 ± 0.4' 


250 ± 154' 








SMM 7 


0.1 ±0.1 


0.3 ± 0.2 


718 ±553 


0.28 ±0.10 


1.5 ± 


0.5 




0.3 ±0.1'' 


0.7±0.3 C 


267 ± 91'' 









" In Paper I, we derived the value iV(N2H + ) « 0.9 1± 0.01 x 10 13 cm 2 from observations of N2H + (1 - 0) emission, which is in good agreement 
with the present value. 
* The principal velocity component was not detected. 
' Calculated by using the optically thin approximation. 
d The associated error is larger than the value. 



yielded a probability of ~ 100%, strengthening the possibility 
that the core masses in these two different parts of Orion B are 
drawn from the same distribution. 



4.4. Virial masses 

In order to study the stability of the cores, we calculated their 
virial masses using the following formula where the effects of 
external pressure and magnetic field are ignored: 



Mw = 



5 fleft'Av^. 

81n2 aG 



(12) 



where G is the gravitational constant, and Av ave is the width of 
the spectral line emitted by the molecule of mean mass ju = 2.33. 
The parameter a = (1 - p/3)/(l - 2p/5), where p is the power- 
law index of the density profile (n(r) oc r~ p ), is a correction for 
deviations from constant density. For starless cores we used the 
value p = 1.0, whereas for protostellar cores we set p = 1.5 (see 
Paper I). The value of M wu calculated with p = 1.0 is about 13% 
higher compared to that calculated with the value p = 1.5. As 
a function of the observed linewidth, Av Q b s (Table |2l Col. (4)), 
Av a ve is given by 



AvL = Av^+Av^t = A V 2 bs + 81n2x 



fcB?kii 

m H 



i --!-i.(i3) 

H yUmol / 



The virial masses are listed in Col. (4) of Table[6] The associated 
error was propagated from those of Av Q b s and 7\j n . 

The virial parameters of the cores were calculated following 
Bertoldi & McKee (1992), i.e., a vir = M vir /M. The uncertainty 
was derived by propagating the errors in both mass estimates. 
The values are given in Col. (5) of Table [6] Note that a vu = 1 
corresponds to the virial equilibrium, 2{T) = -(U), where T 
and U are the internal kinetic and gravitational energies, respec- 
tively. The value a v t = 2 corresponds to the self-gravitating limit 
defined by (T) = -(U). 



5. Discussion 

5.1. The gas kinetic temperature in dense cores in Orion B 

The gas kinetic temperatures within the cores in Orion B9 are in 
the range ~ 9.4 - 13.9 K, with a mean value of 12.0 ± 0.4 K (the 
quoted error is the standard deviation of the mean). The high- 
est temperatures (13.4-13.9 K) are found towards IRAS 05399- 
0121, SMM 4, Ori B9 N, and IRAS 05413-0104. There is, how- 
ever, no such tendency that the warmest cores should be pro- 
tostellar. For example, the Class sources SMM 3 and IRAS 
05405-01 17 have 7\in values at the low end in the sample (11.3 
K for both sources). This is consistent with the notion that em- 
bedded low-mass protostars do not heat significantly the parent 
cores, but the heating is localised in their immediate vicinity 
(< 100 AU; Frie sen et al. 2 010 and references therein). 

Recently, based on the Nobeyama 45-m observations, Ikeda 
et al. (2009) derived the NH3 rotational temperatures towards 
144 positions in Orion B. The cores No. 55-66 in the Ikeda et 
al. (2009) sample lie in the Orion B9 region. For these latter 
cores, they found the values in the range T TO i — 10.1 - 16.1 K. 
By using the relationship ©, this corresponds to the 7\i n values 
10.6 - 18.5 K, with the mean value 16.4 K. This is clearly higher 
than the mean 7\i n of the present study. The difference is likely 
to be caused by the larger beamsize used by Ikeda et al. (~ 1'3). 
We note that the highest value of T\^ n = 43.9 K in the Ikeda et 
al. (2009) sample was found in the NGC 2024 H II region, and 
the lowest values (< 1 1 K) were found in Orion B9 and other 
regions of relatively low star formation activity, i.e., outside the 
NGC 2023, 2024, 2068, and 207 1 regions. 

We note that in several previous submm surveys of 
Orion B dust temperatures around 20 K have been de- 
rived or assumed for dense cores (Johnsto ne et al. 20011 2006; 



Nutter & Ward-Thompson 20071. Together with the results of 



Ikeda et al. (2009) the present temperature determinations 
demonstrate that Orion B contains also several very cold 
cores resembling cores in nearby low-mass star forming re- 
gions (see, e.g., |Rosolowsky et al. 20 08; Sc hnee et al. 2009; 
IFriesen et al. 2009l and references therein). 
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Table 6. Revised core parameters presented in Paper I. 





^eff 


M" 




V 


1' 


N(n 2 y 


<«(H 2 )>" 


Source 


[pc] 


[M~] 


[M Q ] 






[10 22 cm" 2 ] 


[10 4 cm" 3 ] 


IRAS 05399-0121 


0.08 


6.1 ± 1.4 


7.8 ±0.5 


1.3 ± 


0.3 


4.2 ±0.9 


5.5 ± 1.3 


SMM 1 


0.09 


8.4 ± 1.5 


11.1 ±0.3 


1.3 ± 


0.2 


2.8 ± 0.4 


5.3 ±0.9 


SMM2 


0.08 


2.0 ± 0.6 


4.8 ± 0.5 


2.4 ± 


0.8 


1.4 ±0.4 


1.8 ±0.5 


SMM 3 


0.07 


7.8 ± 1.6 


4.0 ± 0.2 


0.5 ± 


0.1 


8.4 ± 1.1 


10.5 ±2.1 


IRAS 05405-0117 


0.07 


2.8 ±0.4 


3.5 ±0.1 


1.2 ± 


0.2 


1.4 ±0.1 


3.8 ±0.5 


SMM 4 


0.07 


2.8 ±0.3 


4.1 ± 0.2 


1.5 ± 


0.2 


1.4 ±0.1 


3.8 ± 0.4 


SMM 4 (2nd v-comp.) 












2.4 ± 0.7 




SMM 5 


0.07 


1.9 ± 0.4 


4.2 ± 0.2 


2.2 ± 


0.5 


1.2 ±0.1 


2.5 ± 0.5 


SMM 6 


0.12 


8.2 ±1.1 


6.2 ±0.2 


0.8 ± 


0.1 


2.0 ±0.1 


2.2 ± 0.3 


Ori B9 N 


0.08 


2.3 ± 0.4 


5.8 ±0.4 


2.5 ± 


0.5 


0.9 ±0.1 


2.1 ± 0.4 


Ori B9 N (2nd v-comp.) 












0.9 ±0.1 




SMM 7 


0.07 


3.6 ± 1.0 


6.9 ± 0.4 


1.9 ± 


0.6 


3.4 ±0.9 


4.8 ± 1.3 


IRAS 05413-0104 


0.06 


2.0 ±0.6 


4.4 ± 0.2 


2.2 ± 


0.7 


3.6 ±0.6 


4.3 ± 1.3 



" The masses and H2 column and number densities were calculated using the values of 7\j„ (= T^t). 
b The table also lists the core virial masses and virial parameters (Cols. (4) and (5)). 



5.2. Kinematics of the core gas 

The NH3 line profiles show that most cores have subsonic non- 
thermal motions (/ tul -b < 1). In two cores (SMM 1 and 7), / tul -b is 
higher than in the rest of the cores, and non-thermal motions ap- 
pear to be slightly transonic (/ tul -b = 1.3 < 2). Recently, Friesen 
et al. (2009) found that the mean / tur b value for dense cores in 
Oph B is 1.5, clearly larger than the corresponding value for our 
cores (a mean / tul -b and its standard deviation is 0.8 + 0.1). 

In Fig. [5] we plot the non-thermal linewidth against the ther- 
mal linewidth for the cores. Figure[5]shows that most of the cores 
in Orion B9 are quiescent, i.e., Avnt < Avj. Also shown in 
this figure are the Avnt - Avj relationships found by Jijina et 
al. (1999) for NH3 cores in clusters and isolated regions. The 
two relationships suggest that non-thermal NH3 linewidths are 
larger for cores in clustered environments than in isolated cores. 
A majority of Orion B9 cores lie on the Jijina plot in the region 
characteristic of isolated cores. 

Like in the case of the kinetic temperature, there is no clear 
difference in linewidths between starless and protostellar cores. 
Protostars are likely to be associated with outflows which could 
cause the linewidths to be larger, but evidence for outflows can 
only be found in IRAS 05399-0121 (HH 92; |Bally et al. 2002) i 
and IRAS 05413-0104 (HH 212; e.g. Xee et al. 2008b . It is pos- 
sible that the HH 92 outflow from IRAS 05399-0121 contribute 
to the relatively broad NH3 lines seen in the adjacent SMM 1 
core which is oriented along the direction of the outflow. 

In Fig. [6] we plot the NH3(1, 1) linewidth as a function of 
core effective radius. Also shown are the linewidth-size relation 
derived by Larson (1981), and the relationship recently found 
by Ikeda et al. (2009) for the H 13 CO + cores in Orion B. The 
Larson relation, thought to arise from the interstellar turbulence, 
was originally presented for the 3D velocity dispersion, cr^,, 
and the maximum linear size of the source, L. In terms of our 
definitions, it can be written as Av[km s -1 ] = 1.9 (R e ff[pc])°' 38 , 
where we have assumed that L = 2R e ff. The linewidths of 
Orion B9 cores mostly deviate downward from the Larson re- 
lation, and the cores lie around the Ikeda et al. (2009) rela- 
tionship, namely Av[kms _1 ] = 1.28 (R core [pc]) ' 38 (we remind 
the reader that the Ikeda et al. sample contains Orion B9 mem- 
bers). On the other hand, it is possible that the more tenuous 
gas in Orion B9 follows the Larson relation: the mean linewidth 
(~ 1.3 km s _1 ) and half-maximum radius (~ 0.6 pc) mea- 



sured from the 13 CO(l - 0) emission towards the Orion B9 re- 
gion by Caselli & Myers (1995) are quite similar to those ex- 
pected from the Larson relation. As pointed out by Maruta et 
al. (2010), this could be explained if the cores were formed in 
regions where supersonic turbulence was dissipated. Dissipation 
of turbulence could be due to shocks in the converging turbu- 
lent flows, where the formation of density enhancements is ex- 
pected to take place (see Sect. 5.9). The timescale for the dis- 
sipation of turbulence in dense cores is comparable to the free- 
fall timescale (e.g., Ma c Low & Klessen 20 04). Thus, because 
the observed linewidth-size relation appears to be quite flat, it 
is possible that the turbulent gas motions at small scales (< 0.1 
pc) are not yet settled into equilibrium state (i.e., relaxed). In 
this case, the core formation should be a rapid process with the 
corresponding timescale being comparable to the free-fall time. 





-0.7 f . 1 ... 1 ... 1 ... 1 ... 1 ... 1 ... ^ 

-0.40-0.38-0.36-0.34-0.32-0.30-0.28-0.26 
og(Av T ) [km s~ ] 

Fig. 5. Non-thermal versus thermal (for H2, i.e., the isothermal 
sound speed) linewidths in log-log scales. The dashed line indi- 
cates Avnt = Avj. The dotted and dash-dotted lines show rela- 
tionships found by Jijina et al. (1999) for NH3 cores with and 
without cluster association, respectively. Symbols have the same 
meaning as in Fig. [4] 
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Fig. 6. Linewidth versus radius for the cores in Orion B9 in log- 
log scales. The Larson (1981) relation is shown by the dotted 
line. The dashed line indicates the relation derived by Ikeda et 
al. (2009) for the H 13 CO + cores in Orion B. Symbols have the 
same meaning as in Fig. [4] 



5.3. Internal pressure support 

The internal kinetic gas pressure, P; nt , within the core consists 
of a thermal pressure, pj = nk^,T^ m , and a non-thermal pressure, 
Pnt = f^mnncrLj, where n is the number density, and fi = 2.33 
(Pint — Pt + Pnt)- In the left panel of Fig. [7] we plot Pi n t/^B 
versus core mass. Excluding the one outlier in the plot (SMM 
6), there is a trend of increasing internal pressure with increas- 
ing core mass. This indicates that the cores are not in pressure 
equilibrium with an external pressure due to, e.g., the weight of 
the cloud in which the cores are embedded dLada et al. 2 008). 
The mean internal pressure of the Orion B9 core population is 
~ 1.1 x 10 6 K cirT 3 . This is in excellent agreement with the 
value ~ 10 6 K cirT 3 estimated by Johnstone et al. (2001) for 
the cores in the northern part of Orion B. For comparison, the 
overall pressure of the ISM in the Galactic midplane, which con- 
sists of thermal and turbulent kinetic pressures, and the pressures 
of magnetic fields and cosmic rays, has been estimated to be 
Pjm/kB ~ 2.8 x 10 4 K cm 3 dBoulares & Cox 19901 1. The high 
pressures inside cores are likely to result from compression by 
gravity (see Sect. 5.4). 

To examine the relative role of turbulence in the core internal 
pressure, we calculated the ratio of thermal to non-thermal pres- 
sure, which is given by R p — 1 1 ff mh (see Sect. 4. 1 .4). The values 
of R p are plotted as a function of mass in the right panel of Fig. [7] 
As can be seen in this figure, thermal pressure is clearly the dom- 
inant source of internal gas pressure for 8 of 1 1 cores. For the 
rest of the cores, thermal pressure support is still significant as 
they have R p > 0.5. Also in this regard, the Orion B9 cores ap- 
pear to be similar to low-mass dense cores in nearby molecular 
clouds, which are commonly found to be thermally dominated 
(e.g., |Myers & Benson 1983llKirk et al. 2007llLada et al. 20081) . 



are self-gravitating (a v i r < 2), and five of them seem to be close 
to the virial equilibrium or collapsing {a wlI < 1). One should 
note that the virial calculation presented above did not include 
any external pressure or magnetic fields. An external presssure 
of Pext/ks ~ 2 - 3 x 10 5 K cirT 3 would bring all starless cores 
except SMM 7 to virial equilibrium or make them collapse (ne- 
glecting the possible magnetic support). Assuming that the gas 
surrounding the cores is characterised by cr NT ~ 0.6 km s _1 and 
n(H2) — 2 x 10 3 cirT 3 , corresponding to the typical width of 
13 CO(l - 0) lines observed by Caselli & Myers (1995) in this re- 
gion and the critical density of the transition in question (which 
is lower than the critical density of the NH3 inversion lines), its 
turbulent ram pressure would roughly equal to the required con- 
fining pressure quoted above. As discussed by Lada et al. (2008) 
in the case of Pipe Nebula, the supersonic intercore turbulence 
can be a manifestation of the self-gravity of the surrounding 
massive cloud. 

It seems likely that most if not all starless cores in Orion B9 
detected in this survey are prestellar, i.e. they will eventually col- 
lapse to stars. The situation resembles that observed in Perseus 
by Foster et al. (2009): most cores are gravitationally bound, and 
~ 1 /3 are in virial equilibrium. 

The virial parameter, a v i r seems to decrease as a function 
of core mass (see Fig. [8] left). The slope of a least-squares fit, 
log(a vir ) = (0.54 + 0.10) -(0.67 + 0.16) log(M), with the linear 
correlation coefficient r = -0.81, is very close to that found by 
Lada et al. (2008) for cores in the Pipe Nebula which are be- 
lieved to be predominantly confined by external pressure. Our 
cores with masses in the range 2-8 M0 correspond to the most 
massive Pipe cores. The slope (~ -2/3) is consistent with the 
theoretical prediction of Bertoldi & McKee (1992) for pressure 
confined cores, and basically results from the fact that there is 
not much core-to-core variation in the average densities and ve- 
locity dispersions. It should be noted, however, that the overall 
level of avir, characterised by the constant of the fit, is lower 
than derived in the Pipe, probably reflecting the different envi- 
ronments and evolutionary stages of the two regions. 

Like Lada et al. (2008) we have also plotted in the right panel 
of Fig. [8]the correlation diagram between the ratio <x 3D /v esc and 
the core mass, where the three-dimensional velocity dispersion, 

cr 3D , is calculated from <x 3D = -^3c;? + 3<x^ T , and the escape ve- 
locity, v esc is given by v esc = y/2GM/R e ff. For all cores, <x 3D is 
smaller than v eS c, supporting the above conclusion that the cores 
are gravitationally bound. Morever, the <x 3D /v esc ratio appears to 
decrease as a function of core mass. A least-squares fit to the 
whole sample gives log (cr 3D /v esc ) = (-0.01 ± 0.03) - (0.43 ± 
0.05) log(M), with r = -0.94. This relationship implies that the 
transition from unbound to bound core occurs at about 0.8-0.9 
M©. This is very close to the peak of the CMF for Orion B 
North, i.e., ~ 1 M© ( |Nutter & Ward-Thompson 2007| l. This is 
expected because the CMFs for Orion B9 and Orion B North 
seem to represent the subsamples of the same parent distribution 
(Sect. 4.3). A similar correspondence was found by Lada et al. 
(2008) for the dense cores in the Pipe Nebula. 



5.4. Dynamical state and gravitational boundedness of the 
cores 

To further examine the dynamical state of the cores, we inspect 
the virial parameters derived in Sect. 4.4. Figure [8] (left panel) 
shows the distribution of the virial parameter for the Orion B9 
cores as a function of mass. Within the errors, all eleven cores 



5.5. Correlating the properties of protostellar cores 

In Paper I, we derived the SEDs for the protostellar cores in 
Orion B9. For all sources, the observed flux densities were fit- 
ted by a two-temperature (warm+cold) composite model. The 
bolometric dust temperatures and bolometric luminosities of the 
sources derived from the SEDs are given in Table 6 of Paper I. 
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Fig. 7. Left: Total internal kinetic gas pressure versus core mass. Right: The ratio of thermal to non-thermal pressure as a function 
of mass. The dashed horizontal line marks the point at which thermal and non-thermal pressure support are equal. Both plots are in 
log-log scales, and the symbols have the same meaning as in Fig. |4] 




Fig. 8. Left: Virial parameter versus core mass. The dashed line indicates the virial equilibrium limit (a vll - = 1), and the dash-dotted 
line shows the limit of gravitational boundedness (a v i r = 2). Right: Ratio of 3D velocity dispersion to escape velocity against core 
mass. The dashed line indicates cr3 D /v eS c = 1. All cores are likely to be gravitationally bound. The solid line in both plots indicates 
the least-squares fit to the data. Both plots are in log-log scales, and the symbol key is identical to that of Fig. [4] 



In the left panel of Fig. [9] we show the correlation plot be- 
tween Tidn and the bolometric temperature, 7b i- Only for SMM 
3, the two values are similar within the errors, and in the case of 
SMM 4, 7\j n and T\, \ agree within a few kelvins. For the other 
sources (i.e., the three IRAS sources), however, the temperatures 
estimated from the SEDs are clearly higher than those obtained 
from NH3 measurements. The SEDs of SMM 3 and 4 were con- 
structed by using only three flux density values at 24, 70, and 
870 jt/m (see Paper I; Fig. 6 therein). In the case of other sources, 
also the flux densities at all four IRAS bands (12, 25, 60, and 
100 fj.m) were included in the SED fits. This may have caused 
the temperature of the cold part of the spectrum to be overesti- 
mated. Shetty et al. (2009) recently concluded that if using flux 
density values at widely separated wavelengths, and including 
short-wavelength flux densities in the SED fit, the obtained tem- 
peratures may be too high. 

The right panel of Fig. [9] shows a plot of Avnt as a func- 
tion of Lboi- There appears to be a positive correlation between 
the two quantities, and the least-squares linear fit to the data 
gives log(Av NT ) = (-0.65 ± 0.10) + (0.25 ± 0.11) log(L bo i), with 



r - 0.80. We also show the relationships found by Jijina et al. 
(1999) for NH3 cores in clusters and cores without cluster asso- 
ciation (see their Table B7), and the Myers et al. (1991) relation 
for a sample of 61 IRAS point sources (their Eq. (3)). Our re- 
sult is similar to the relationship found in the latter two studies, 
Avnt k ^boi 9 - The fact that Avnt increases as a function of Lboi 
can be explained by the effects of winds driven by the embedded 
protostar+disk system, and by the fact that more turbulent initial 
conditions give rise to more massive, i.e., more luminous stars 
(Myers et al. 1991 1. As in the case of Avnt - Avt relation (see 
Fig. [5}, the Orion B9 cores more closely follow the relationship 
for isolated cores. 



5.6. Column densities and fractional abundances of A/H 3 
N 2 H + 



and 



The mean NH3 fractional abundance in our sample is 4.8 + 0.7 X 
10" , where the ±-error represents the standard deviation of the 
mean. The value of (jt(NH3)) in Orion B9 is similar to those 
recently found by Foster et al. (2009) and Friesen et al. (2009) 
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Fig. 9. Left: T^ n versus T\, \. The dashed line indicates 7\j n = 7boi- Right: Avnt versus Lboi- The solid line indicates the least- 
squares fit to the data. The upper and lower dashed lines represent relationships found by Jijina et al. (1999) for NH3 cores with and 
without cluster association, respectively. The dash-dotted line shows the relationships found by Myers et al. (1991) for IRAS point 
sources. Both plots are in log-log scales. 



for the dense cores in Perseus and Ophiuchus, respectively (i.e., 
a few to several xlCT 8 ). For comparison, the NH3 abundances 
in dense cores in nearby regions of isolated star formation, such 
as Taurus, are commonly found to be ~ 0.4 - 5 x 10~ 8 (e.g., 
IHotzel et al. 200TllTafalla et al. 20021 2006). The average N 2 H + 
abundance in Orion B9, <*(N 2 H + )> = 2.5 ± 0.8 x 10" 10 , is very 
similar to those in isolated low-mass cores (e.g., (.x(N 2 H + )) — 
2+ 1 x 10~ 1() and 3+2x 10~ 10 for the starless and protostellar core 
samples of Caselli et al. (2002)). On the other hand, the value of 
(jc(N 2 H + )) is about 2-6 times lower than those determined by 
Friesen et al. (2010) in the Oph Bl and B2 cores which lie in the 
region of clustered star formation. 

As shown in the left panel of Fig. [10] the fractional NH3 
abundance appears to decrease with increasing H2 column den- 
sity (as traced by 870 /mi continuum emission smoothed to 
40" resolution). A linear-least-squares fit to the data gives 
log *(NH 3 ) = (4.2 + 6.6) - (0.5 ± 0.3) logjV( H 2 ), wit h r = -0.5 1 . 
The slope of this relationship, ^(NHj) <x 1/ VAf(H 2 ), is very sim- 
ilar to that found by Friesen et al. (2009, 2010) in the Oph B2 
core. In the middle panel of Fig. [10] we plot x(NH3) as a func- 
tion of «(H 2 ). The trend is similar as above, and the least-squares 
lineai- fit to the data indicates log jc(NH 3 ) = (-4.5 ± 1.4) - (0.6 + 
0.3) logn(H 2 ), with r = -0.57. The tendencies described above 
are likely to manifest the fact that NH3, like species contain- 
ing carbon or oxygen, accrete onto grain surfaces at high den- 
sities. According to chemistry models of Aikawa et al. (2005) 
and Flower et al. (2006), the NH3 abundance becomes heavily 
depleted at densities «h Z a few xlO 6 crrT 3 (due to freeze-out of 
the parent species N 2 ). The x(NH3) - n(H 2 ) relationship shown 
in Fig. [10] agrees reasonably well with the predictions for the 
early stages of collapse presented in Flower et al. (2006; see, 
e.g., their Fig. 4). In the Flower et al. model, the protostellar col- 
lapse occurs on the free-fall timescale which is about 4 x 10 5 yr 
at «h = 10 4 cm 4 . This is in agreement with statistical estimates 
of the prestellar core lifetime (see Sect. 5.9). 

Unlike in Friesen et al. (2010), no correlation was found be- 
tween x(N 2 H + ) and H 2 column or number densities. In the right 
panel of Fig.. QJJ we plot the NH3/N 2 H + column density ratio as 
a function of A^(H 2 ) (determined using 40" resolution). A pos- 
sible negative correlation can be seen between the two quanti- 
ties, resembling the result found by Friesen et al. (2010) in Oph 



Bl, and by Johnstone et al. (2010) in Perseus. As can be seen 
from the figure and Col. (4) of Table[5] prestellar cores generally 
have higher values of A^(NH3)/A^(N 2 H + ) ratio than in protostel- 
lar cores. The tendency has been observed also in other star- 
forming regions (Paper I and references therein; Friesen et al. 
2010). The reason for this differentiation is not quite clear. The 
abundances of both molecules build up slowly, and they benefit 
from increasing density and the freezing of CO and other heav- 
ier species (e.g. JAikawa et al. 20051) . For both molecules the en- 
hanced formation at high densities should be counteracted by 
accretion onto dust grains. In the time evolution of the column 
densities these contradictory tendencies are probably reflected 
as NH3 and N 2 H + peaks coming one after the other. Because of 
the different radial distributions of these molecules (and proba- 
bly also because of the different critical densities of their tran- 
sitions) the derived column density ratios depend on the spatial 
resolution. The fact that Friesen et al. (2010) could see the effect 
of N 2 H + depletion towards the highest densities can probably 
be explained by a higher resolution and better data quality than 
available here. 

Our data contains, however, evidence that N 2 H + is frozen 
out in the core SMM 4, as suggested in Paper I. The core was 
not detected in N 2 H + (3 - 2), and the absence of this molecule is 
most probably not caused by efficient desorption of CO, which 
is one of the principal destructors of N 2 H + in the gas phase. 
The gas kinetic temperature derived from ammonia (~ 14 K) 
is clearly lower than the sublimation temperature of CO (~ 20 
K; I Aikawa et al. 20081 1. Depletion of N 2 H + could also be the 
reason for the very weak N 2 H + (3 - 2) emission seen towards 
Ori B9 N. An other indication of depletion in the core Ori B9 
N is that the H 2 D + (li - In) line at ~ 9 km s 1 is detected 
in this source, about 38"southeast from the dust peak position 
( Harju et al. 2006] ). This molecule is known to resist depletion 
"to last" and is observed in highly depleted cold cloud cores 
dCaselli et al. 20081 

5. 7. Low-velocity gas emission seen towards Orion B9 

The LSR velocities of SMM 7 (3.6 km s 1 ) and IRAS 05413- 
0104 (1.5 km s _1 ) are about 6-8 km s _1 less than the systemic 
velocity of Orion B (see Sect. 3.1). We note that the systemic 
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Fig. 10. Comparison of chemical abundances as a function of H2 column and number density. All plots are in log-log scales, and 
the symbol key is identical to that of Fig. [4] The solid lines indicate the least-squares fits to the data. Left: x(NH3) versus N(H2). 
The dashed line indicates the relationship found by Friesen et al. (2009, 2010) for Oph B2. Middle: x(NH 3 ) versus n(H 2 ). Right: 
NH3/N2I-F column density ratio versus N(H2). 



velocity of IRAS 05413-0104 is known to be low before the 
present study (e.g., from NH3 measurements by Wiseman et al. 
(2001)). Given that the sources north-east from the central part 
of Orion B9 appear to have lower radial velocities, the same is 
likely to be true for IRAS 05412-0105. Indeed, IRAS 05412- 
0105 and 05413-0104 could be embedded in a common 1.2 mm 
clump (Kauff mann et al. 2008i l. No NH 3 (1, 1) line was detected 
towards IRAS 05412-0105 (Sect. 2.1), and we could not found 
any other spectral line observations for it from the literature (pro- 
viding kinematic information). Also, the second velocity compo- 
nents seen in the NH3 and N2H + (3 - 2) spectra towards SMM 
4 and Ori B9 N have centroid velocities ~ 7 - 8 km s _1 lower 
than the systemic velocity. In Paper I, the ~ 2 km s _1 velocity 
components were detected in N2H + (1 - 0) and/or N2D + (2 - 1) 
towards a few target positions near IRAS 05405-0117 (called 
Ori B9 E in Paper I) and Ori B9 N. Similarly, Aoyama et al. 
(2001) found that the C ls O(l - 0) and H 13 CO + (l - 0) clumps 
associated with IRAS 05405-01 17 have centroid LSR velocities 
of about 2 km s _1 (at the resolutions 2.'7 and 3'8, respectively). 
This raises the question whether the sources and/or high-density 
gas showing the lower velocity emission are associated with the 
same cloud complex as the "regular" 9 km s _1 component ? 

The CO(l - 0) integrated intensity maps in the velocity 
ranges -1.0 - 2.3 and 2.3 - 5.5 km s" 1 by Wilson et al. (2005; 
their Fig. 3) show relatively strong emission in the direction of 
Orion B9. As discussed by Wilson et al. (2005), the stellar winds 
from the Ori OB lb subgroup have likely interacted with the low- 
longitude and -latitude part of Orion B, resulting in an acceler- 
ated motion towards the Sun. This fraction of the gas is likely to 
be located a few tens of parsecs closer to the Sun than the higher 
velocity gas, and it is estimated to have a mass of about 15% of 
the total mass of Orion B, i.e., ~ 10 4 M Q dWilson et al. 2 005). 
Our pointed molecular line observations in high density tracers, 
including deuterated species, show that the low-velocity gas in 
Orion B also contains dense cores, possibly originating from the 
fragmentation of the cloud region compressed by the feedback 
from the massive stars of the Ori OB lb group. 

5.8. Core spatial distribution revisited 

In Paper I, we studied the spatial distribution of cores in Orion 
B9 in order to examine the fragmentation length-scale. We de- 
termined the core-separation distribution and the number dis- 
tribution of the projected separation distance between nearest 
neighbours. These were compared with the corresponding ran- 



dom distributions. We concluded that the observed distributions 
are random-like, and thus the origin of cores is possibly caused 
by turbulent fragmentation; random distribution is expected if 
the cloud fragmentation is driven by a stochastic turbulent pro- 
cess. As discussed in the previous subsection, some of the cores 
are probably located at somewhat different distance than the "9 
km s" 1 members of Orion B9". For this reason, we re-examined 
the core spatial distribution in the region. 

We note furthermore that there was a flaw in our IDL proce- 
dure used to generate the histograms shown in Fig. 8 of Paper I. 
The number of separations (66 for the core separations, and 12 
for the nearest neighbours) were not correct, and the numbers in 
the observed and model (random) distributions were not equal. 
Here we reproduce the histograms and show the new plots for 
the core separations and nearest neighbours in the top panels of 
Fig-El In the bottom panels of Fig.[TT] we show the same dis- 
tributions but without the sources SMM 7, IRAS 05412-0105, 
and IRAS 05413-0104. In the latter figure, the size of the region 
used to generate the random distributions was decreased from 
0.22 n° to 0.06 D°. The obtained statistics of the distributions 
are presented in Table [7] In this table, we give the mean (and 
its standard deviation) and median of the observed core spatial 
distribution (Cols. (2) and (3)), those of the corresponding ran- 
dom distribution (Cols. (4) and (5)), ratios between the observed 
and random mean and median separations (Cols. (6) and (7)), 
and probability given by the two-sample K-S test that the ob- 
served and random distributions are drawn from the same un- 
derlying distribution, see below (Col. (8)). The random distri- 
butions were generated a hundred times and the averaged his- 
tograms were compared with those derived from observations. 
The numerical values for the random distributions presented in 
Table [7] are the averages and the standard deviations of these 
100 runs. Only a few values differ from those reported in Paper 
I (see the footnotes in the table). Also, for the core positions 
in Orion B North (Nutter & Ward-Thompson 2007), the values 
are the same as those reported in Paper I. We note that for the 
core-separation distributions (r)o r iB9/('")oriBN = 0.63 ± 0.01 and 
'briB9/'briBN = 0.62, and for the nearest-neighbour distribution 
<r>OriB9/<>">OriBN = 1-95 ± 0.27 and r riB9/roriBN = 1-87. 

In Paper I, we suggested that the comparable mean and me- 
dian values between the observed and random distributions in- 
dicate that the cores are likely to be randomly distributed within 
the region. To examine this in more detail, we carried out K-S 
tests between the observed and model distributions. As shown 
in the last column of Table [7] the probability that the observed 
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distribution and the generated random distribution represent the 
same underlying distribution is very small in all cases. Such low 
probabilities call into question the similarity between the ob- 
served and random distributions. Thus, even the observed and 
model mean and median separation-distances are comparable, 
the above K-S test probabilities suggest that the distributions as 
a whole are not similar. The K-S probability is expected to be 
a more robust measure of the similarity between the two dis- 
tributions than the comparison of the mean and median values 
because the latter two can be the same for two different distribu- 
tions. 

The observed mean and median distances between the 
nearest-neighbours are in the range ~ 3 x 10 4 - 6x 10 4 AU. 
Assuming that the parental cloud region is characterised by the 
average gas kinetic temperature and density values of 15 K and 
~ 2 x 10 3 - 10 4 cm -3 , respectively (Sect. 5.4), the thermal Jeans 

length, Aj - yfircf/Gp, is about 5 x 10 4 - 10 5 AU. These are com- 
parable to the observed core separations indicated above. This 
suggests that the fragmentation of the region into cores is caused 
by gravitational instability. The parental cloud region in which 
the cores have formed could have been initially compressed by 
the winds from the nearby massive stars as discussed in Sect. 5.7. 



the dominance of turbulent fragmentation in core formation. 
Moreover, we stated that because region represents clustered star 
formation, the turbulence is most likely driven on large scales 
(e.g., IKlessen et al. 20001 IKlessen 20011 >. However, the present 
molecular line data and a more careful analysis of the spatial dis- 
tribution have called both arguments into question: According to 
the K-S test, the core distribution is unlikely to be random, and 
the core properties resemble closely those in the regions of iso- 
lated star formation. 

To summarise, the current data suggest that the evolution of 
the Orion B9 cores is dynamic. As discussed in Sect. 5.7, the 
Orion B9 region is likely to be influenced by the massive stars 
of the Ori OB lb group. This interaction process could have led 
to the driving of the large-scale turbulence and/or to the com- 
pression of the cloud material into dense sheets and filaments as 
observed in molecular clouds. These structures could have then 
been fragmented into dense cores via gravitational instability as 
suggested by the Jeans-length analysis in Sect. 5.8. Further in- 
vestigations of the chemistry, deuteration, and the degree of ion- 
isation can help to fix the evolutionary timescales, and thereby 
estimate the roles of turbulence, thermo-gravitational instability, 
and ambipolar diffusion in the core formation. 



5.9. The lifetime and origin of dense cores in Orion B9 



6. Summary and conclusions 



According to the scenario of turbulence-regulated star for- 
mation, the prestellar core evolution is dynamic and the 
corresponding lifetime is only a few times the free-fall time 
dMac Low & Klessen 20041 |Vazquez-Semadeni et al. 2005 ). 



Several results of our study conform to this scenario: as was 
discussed in Sect. 5.2, the Orion B9 cores could have been 
formed in turbulent shocks where the dissipation of the kinetic 
energy of turbulent motions took place. It was also concluded 
that the timescale for core formation should be comparable to 
the free-fall time. On the other hand, the abundances of the 
N-bearing species studied in the present paper are consistent 
with those predicted by chemical models based on dynamical 
prestellar core evolution (Sect. 5.6). Moreover, in Paper I we 
deduced that there are equal numbers of pre- and protostellar 
cores in the Orion B9 region, which suggest that the corre- 
sponding lifetimes are comparable. Based on similar statistical 
results, i.e., N ps:e /Np Ioto ~ 1, recent studies have shown that the 
duration of the prestellar phase of core evolution at densities 
nQiz) Z 10 4 cirT 3 is only a few times the free-fall time (see 
references in Paper I; Ev ans et al. 2 009). It is possible that some 
of the cores do not belong to the same volume of the cloud 
as those with centroid velocities around 9 km s _1 . Even if the 
cores with lower centroid velocities are excluded (i.e., SMM 
7, and IRAS 05413-0104 and 05412-0105), the above ratio is 
5:4= 1.25, which is still consistent with the dynamic core 
evolution scenario. 

In the present study (and in Paper I), we found that the 
CMFs in Orion B9 and Orion B North very likely represent the 
subsamples of the same parent distribution, and resemble the 
stellar IMF ( |Nutter & Ward-Thompson 2QQ7\ . Numerical simu- 
lations of turbulent fragmentation have been able to reproduce 
the general shape of the IMF (e.g., IPadoan & Nordlund 20021 
Ballesteros-Paredes et al. 20061 1. On the other hand, recent mod- 
els concerning fragmentation initiated by the ambipolar diffu- 
sion (AD) have likewise succeeded in explaining the connection 
between the CMF and the IMF (IKunz & Mouschovias 20091 1. 

In Paper I, we used the spatial distribution of cores, which 
seemed to mimic a random distribution, as an argument for 



We observed the NH.3(1, 1) and (2, 2) inversion transitions, and 
the N2H + (3 - 2) lines towards the submm cores in Orion B9 
star-forming region with the Effelsberg 100-m and APEX tele- 
scopes, respectively. These line observations were combined 
with our previous 870 //m submm dust continuum data of the 
region. The data were used to determine the physical char- 
acteristics of dense cores in this region. We mainly investi- 
gated the temperatures and kinematics of the cores, and recal- 
culated the temperature-dependent quantities presented in Paper 
I (Miettin en et al. 20 09). Our main results and conclusions are 
as follows: 

1 . The gas kinetic temperatures of the cores are in the range ~ 
9.4- 13.9 K, with an average value of about 12 K. No significant 
difference in T^ n was found between starless and protostellar 
cores. The temperature values are similar to those in nearby low- 
mass star-forming regions. 

2. The cores are characterised by subsonic, or at most tran- 
sonic, non-thermal motions. Thus, the cores are kinematically 
similar to those in nearby low-mass star-forming regions. 

3. In the case of protostellar cores, we found a positive cor- 
relation between the non-thermal NH3 linewidth and bolometric 
luminosity. This can be understood so that on one hand, the em- 
bedded central protostar enhances the level of turbulent motions 
within the parent core, and on the other hand, more turbulent ini- 
tial conditions lead to the formation of more massive stars with 
higher luminosity. 

4. The core masses are, on average, about 11% higher than 
those reported in Paper I when calculated by assuming that 
?dust = Tkin (in Paper I, we assumed that Tjust = 10 K for starless 
cores). These masses are very likely drawn from the same parent 
distribution as the core masses in the northern part of Orion B 



UNutter & Ward-Thompson 2007) . 

5. Almost all cores were found to be close to virial equi- 
librium when the internal kinetic pressure and self-gravity are 
taken into account. In addition, they are likely to be exposed to 
a substantial turbulent ram pressure from the intercore gas (as 
can be evaluated from previous 13 CO data), which is sufficient 
to bring them to equilibrium or make them collapse. In particu- 
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Fig. 11. Top left: Observed core-separation distribution (solid line) compared with the expected distribution for random distribution 
of the same number of cores as the observed sample over an identical area (dashed line). Top right: Same as in the top left panel 
but for the nearest-neighbour distribution. Bottom panels: Same as in the top panels, but SMM 7, IRAS 05412-0105, and IRAS 
05413-0104 have been excluded from the samples. 



Table 7. Statistics of the core spatial distributions in Orion B9. 
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" In Paper I, this was reported to be 5.59 ± 0.05. 

* SMM 7, IRAS 05412-0105, and IRAS 05413-0104 have been excluded from the sample. 

c 5.08 ± 0.08 in Paper I. 

d 5.07 ±0.11 in Paper I. 



lar, it seems probable that most if not all starless cores detected 
in this survey are prestellar. 

6. The fractional NH3 and N2H + abundances in the cores 
are ~ 1.5 - 9.8 x 10~ 8 and ~ 0.2 - 5.9 x 10~ 10 , respectively. 
The NH3/N2H + column density ratio is higher in prestellar cores 
than protostellar cores. Similar trend has also been observed in 
other star-forming regions. The NH3 abundance appears to de- 
crease with increasing gas density (as traced by 870 /im contin- 
uum emission). This tendency is likely to be caused by the ac- 



cretion of NH3 onto grains at very high densities, and conforms 
to the recent results of Friesen et al. (2009) in Ophiuchus. 

7. A few cores have a much lower radial velocity compared 
to the systemic velocity of the region (~ 9 km s _1 ), and do not 
necessarily belong to the same volume of the cloud as the rest 
of the cores. Instead, they are likely to be members of the "low- 
velocity" part of Orion B. This increases the relative number of 
pre- and protostellar cores in Orion B9 and consequently, the 
statistical lifetime of the prestellar phase of core evolution (by 
25% in both cases). 
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8. Many of the properties of dense cores in Orion B9 suggest 
that their evolution is likely to be dynamic, i.e., comparable to 
the free-fall timescale. On the other hand, the new data and an 
improved analysis of spatial distribution have made the role of 
turbulent fragmentation in the core formation less evident than 
it appeared in Paper I. According to the new results, the cores 
are unlikely to be randomly positioned within the region, and 
they are not clustered in the way expected if the fragmentation 
is driven on large scale turbulence. Futhermore, the core prop- 
erties are similar to those found in the regions of isolated star 
formation. 

9. The Orion B9 region is likely to be influenced by the mas- 
sive stars of the nearby Ori OB lb group. This explains the ori- 
gin of the low-velocity material in the southern parts of Orion 
B dWilson et al. 2005I) . Moreover, this interaction process could 
have led to formation of dense filaments, from which the dense 
cores were fragmented out by the action of gravitational insta- 
bility. 

The present paper demonstrates the importance of a kine- 
matic and temperature information for the studies of dense cores 
and star formation. A lack of radial velocity information, for in- 
stance, may cause errors when examining the relative numbers 
of starless and protostellar cores. This, in turn, affects the life- 
time estimates based on statistical arguments. The core lifetime 
is an important discriminator between different theories of core 
and star formation, particularly between the turbulence and am- 
bipolar diffusion driven star formation scenarios. 
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